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We discuss the production of charmonium and open charm in nuclear colhsions at SPS/FAIR 
energies within the framework of the statistical hadronization model. The increasing importance 
at lower energies of production is discussed and provides a challenge for future experiments. 
We also show that possible modifications of charmed hadrons in the hot hadronic medium do not 
lead to measurable changes in the cross sections for D-meson production. A possible influence 
of medium effects can be seen, however, in yields of charmonium. These effects are visible at all 
energies and results will be presented for the energy range between charm threshold and RHIC 
energy. 
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1. Introduction 



Charmonium production is considered, since the original proposal more than 20 years ago 
about its suppression in a Quark-Gluon Plasma (QGP) [|l|], as an important probe to determine 
the degree of deconfinement reached in the fireball produced in ultra-relativistic nucleus-nucleus 
collisions. In the original scenario of J/i//^ suppression via Debye screening [[T]] it is assumed that 
the charmonia are rapidly formed in initial hard collisions but are subsequently destroyed in the 
QGP (see an update of this picture in ref. [|2|]). 

In a recent series of publications ^ ^, ^] we have demonstrated that, in the energy range 
from top SPS energy {^/snn ~ 17 GeV) on, the data on J/y and y' production in nucleus-nucleus 
collisions can be well described within the statistical hadronization model proposed in This 
includes the centrality and rapidity dependence of recent data at RHIC {^ysNN=200 GeV) published 
by the PHENIX collaboration [Q]. We note that the extrapolation of these results to LHC energy 
i\/^NN=5-5 TeV) yields a rather striking centrality dependence ||^ ^. Depending on the magnitude 
of the cc cross section in central Pb-Pb collisions even an enhancement of J/i/a production 
compared to pp collisions (R^/^ > 1) is expected due to hadronization (at chemical freeze-out) of 
uncorrected (at these high energies) charm quarks thermalized in QGP 

Here we explore the lower energy range (from near threshold, ^^smn ~ 6 GeV), which can 
be investigated in the CBM experiment [||] at the future FAIR facility. One of the motivations for 
such studies was the expectation [T^ to provide, by a measurement of D-meson production near 
threshold, information on their possible in-medium modification near the phase boundary. How- 
ever, the cross section Ocs is governed by the mass of the charm quark /Mc ~ 1.3 GeV, which is much 
larger than any soft Quantum Chromodynamics (QCD) scale such as Aqcd- Therefore we expect 
no medium effects on this quantity. ^ The much later formed D-mesons, or other charmed hadrons, 
may well change their mass in the hot medium. The results of various studies on in-medium 
modification of charmed hadrons masses [IC, 11, 12, [l^, [l^, |l^, L^l sometimes con- 
tradictory. Whatever the medium effects may be, they can, because of the charm conservation, 

Ccc = f (cd + <^A, + <yz, + ■■■) + {<^ri, + <^j/\if + (^Xc ~^ ■■■)' ^'^ly 1^^^ ^'^ ^ redistribution of charm 
quarks [19]. This argument is essentially model-independent and applies equally well at all ener- 
gies. Here we will consider various types of scenarios for medium modifications and study their 
effect within the statistical hadronization framework in the energy range from charm threshold to 
collider energies. In this context, we note that excellent fits of the common (non-charmed) hadrons 
to predictions of the thermal model have been obtained using vacuum masses (see ref. [ ^0| ] and 



references therein). An attempt to use modified masses for the RHIC energy |21] has not produced 
a conclusive preference for any mass or width modifications of hadrons in medium. On the other 
hand, some evidence for possible mass modifications was presented in the chiral model of [22]. 



2. Assumptions and ingredients of the statistical hadronization model 

The statistical hadronization model (SHM) [^, ^ assumes that the charm quarks are produced 
in primary hard collisions and that their total number stays constant until hadronization. Another 

' Such a separation of scales is not possible for strangeness production, and the situation there is not easily compa- 
rable. 
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important factor is thermal equilibration in the QGP, at least near the critical temperature, Tc. We 
neglect charmonium production in the nuclear corona [Q], since we focus in the following on central 
collisions {Npart=i'^^), where such effects are small. 

In the following we briefly outline the calculation steps in our model [||, ^ . The model has the 
following input parameters: i) charm production cross section in pp collisions; ii) characteristics 
at chemical freeze-out: temperature, T , baryochemical potential, ptb, and volume corresponding to 
one unit of rapidity V^y=\ (our calculations are for midrapidity). Since, in the end, our main results 
will be ratios of hadrons with charm quarks nomalized to the cc yield, the detailed magnitude of the 
open charm cross section and whether to use integrated yield or midrapidity yields is not crucial. 
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Figure 1: Energy dependence of the 
charm production cross section in pp 
collisions. The NLO pQCD values 
[ p4| ] are compared to calculations us- 
ing PYTHIA and to data in pA col- 
lisions, taken from ref. [|5]|. Our 
extrapolations for low energies are 
shown with continuous lines, for to- 
tal and midrapidity (ddtr/dy) cross 
section. The open square is a midra- 
pidity measurement in pp collisions 
p^ ] . The dashed line with dots indi- 
cates a parameterization of the mea- 
sured energy dependence of the 
production cross section [EtI]. 



The charm balance equation [g], which has to include canonical suppression factors [23] when- 
ever the number of charm pairs is not much larger than I, is used to determine a fugacity factor gc 
via: 



n: 



dir 



2^' "'kigcN'J^c) 



■th 



(2.1) 



Here N^l^^ is the number of initially produced cc pairs and !„ are modified Bessel functions. In 



the fireball of volume V the total number of open {N'^l = ri^^W) and hidden (A^^^ = r{^^/) charm 
hadrons is computed from their grand-canonical densities r^^^ and n%, respectively. This charm 
balance equation is the implementation within our model of the charm conservation constraint. The 
densities of different particle species in the grand canonical ensemble are calculated following the 
statistical model [20|. The balance equation (2.1) defines the fugacity parameter g^- that accounts 
for deviations of heavy quark multiplicity from the value that is expected in complete chemical 
equilibrium. The yield of charmonia of type j is obtained as: N j = g^N'-^, while the yield of open 
charm hadrons is: A^,- = gcNj^h{gcK'',)/Io{gcKc)- 
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Vs^N (GeV) Vs^^ (GeV) 

Figure 2: Energy dependence of the number of ini- Figure 3: Energy dependence of the volume at 
tially produced charm quark paks (Npart=^50). midrapidity, Vav=1i for central collisions [^. 



As no information on the charm production cross section is available for energies below v^=15 
GeV, we have to rely on extrapolation. The basis for this extrapolation is the energy dependence of 
the total charm production cross section calculated in ref. [ p4| ] for the CTEQ5M parton distribution 
functions in next-to-leading order (NLO), as shown in Fig. [l|. We have scaled these calculations to 
match the more recent values calculated at ^/s=200 GeV in ref. [28]. We employ a threshold-based 
extrapolation using the following expression: 

Occ = kil-^r/VsnV^-/V~sf (2.2) 



with A:=1.85 /^b, y^5^=4.5 GeV (calculated assuming a charm quark mass mc=l3 GeV []29[]), 
a=4.3, and b=-l.44. The parameters a, b, k were tuned to reproduce the low-energy part of the 
(scaled) NLO curve. The extrapolated curves for charm production cross section are shown with 
continuous lines in Fig. ^. Also shown for comparison are calculations with PYTHIA [25|. To ob- 
tain the values at midrapidity we have extrapolated to lower energies the rapidity widths (FWHM) 
of the charm cross section known to be about 4 units at RHIC [ ^ ] and about 2 units at SPS [pO|]. 
With these cross section values, the rapidity density of initially produced charm quark pairs, shown 
in Fig. m strongly rises from 1.1-10"^ to 1.7 for the energy range ^^^^=7-200 GeV. We note that 
the so-obtained charm production cross section has an energy dependence similar to that measured 



for production, recently compiled and parametrized by the HERA-B collaboration [g7|]. For 
comparison, this is also shown in Fig. [l[ The extrapolation procedure for the low-energy part of 
the cross section obviously implies significant uncertainties. We emphasize, however, that the most 
robust predictions of our model, i.e. the yields of charmed hadrons and charmonia relative to the 
initially produced cc pair yield are not influenced by the details of this extrapolation. 

For the studied energy range, ^5^^=7-200 GeV, T rises from 151 to 161 MeV from ^snn=7 
to 12 GeV and stays constant for higher energies, while Hh decreases from 434 to 22 MeV [EH]. 
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Vs^N (GeV) Vs^^ (GeV) 

Figure 4: Energy dependence of the canonical suppression for charm, I\ /Iq (left panel) and of the charm 
quark fugacity, gc (right panel). 



The volume Vaj,=i at midrapidity, shown in Fig. || [ pO| ] continuously rises from 760 to 2400 fm^. 
Due to the strong energy dependence of charm production, Fig. ^ the canonical suppression factor 
ih/Io) varies from 1/30 to 1/1.2. Correspondingly, the charm fugacity gc increases from 0.96 to 
8.9, see Fig. |. 

Before proceeding to discuss our results, we would like to emphasize some peculiar aspects 
of charm at low energies. First, the assumption of charm equilibration can be questionable. In this 
exploratory study we have nevertheless assumed full thermalization. At SPS and lower energies 
collision time, plasma formation time, and charmonium (or open charm hadrons) formation time 
are all of the same order |31]. Furthermore, the maximum plasma temperature may not exceed the 
J/Y dissociation temperature. To, although recent results [32] indicate that To can be very close 
to Tf. Charmonia may be broken up by by gluons and by high energy nucleons still passing by 
from the collision. In this latter case cold nuclear suppression needs to be carefully considered (as 
discussed, e.g., in [33, 34]). Consequently, our calculations, in which both charmonium formation 
before QGP production and cold nuclear suppression are neglected, may somewhat underestimate 
the charmonium production yield at SPS energies [Q] and below. 

We note that models that combine the 'melting scenario' with statistical hadronization have 



been proposed [|35|]. Alternatively, charmonium formation by coalescence in the plasma [ |36| , |37 
38, ^ as well as within transport model approaches [p3, 41] has been considered. 



3. Energy dependence of charmed hadrons yield 

Our main results are presented in Fig. ^ The left panel shows our predictions for the energy 
dependence of midrapidity yields for various charmed hadrons. Beyond the generally decreasing 
trend towards low energies for all yields one notices first a striking behavior of the production of A+ 
baryons: their yield exhibits a weaker energy dependence than observed for other charmed hadrons. 
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In our approach this is caused by the increase in baryochemical potential towards lower energies 
(coupled with the charm neutrality condition). A similar behavior is seen for the S+ baryon. These 
results emphasize the importance of measuring, in addition to D-meson production, also the yield 
of charmed baryons to get a good measure of the total charm production cross section. In detail, 
the production yields of D-mesons depend also on their quark content. 




Figure 5: Energy dependence of charmed hadron production at midrapidity. Left panel: absolute yields, 
right panel: yields relative to the number of cc pairs. Note, in both panels, the scale factors of 10 and 100 
for y/i/ and y/' mesons, respectively. 

The differing energy dependences of the yields of charmed hadrons are even more evident in 
the right panel of Fig. ^, where we show the predicted yields normalized to the number of initially 
produced cc pairs. Except very near threshold, the //i/A production yield per cc pair exhibits a slow 
increase with increasing energy. This increase is a consequence of the quadratic term in the J/i/a 
yield equation discussed above. At LHC energy, the yield ratio J/xjf/cc approaches 1% [HI], scaling 



linearly with Occ (for details see [|T9|]). The y' yield shows a similar energy dependence as the 
except for our lowest energies, where the difference is due to the decrease of temperature 
(see above). We emphasize again that this model prediction, namely yields relative to cc pairs, is a 
robust result, as it is in the first order independent on the charm production cross section. Due to the 
expected similar temperature, the relative abundance of open charm hadrons at LHC is predicted 
[^] to be similar to that at RHIC energies. 

4. Effects of in-medium modification of charmed hadrons masses 

We consider two scenarios^ for a possible mass change Am of open charm hadrons containing 
light, u or d, quarks: i) a common decrease of 50 MeV for all charmed mesons and their antiparti- 

^The scenarios are constructed by modification of the constituent quark masses of light (u and d) quarks in the 
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cles and a decrease of 100 MeV for the Ac and Zc baryons (50 MeV decrease for S^); ii) a decrease 
of 100 MeV for all charmed mesons and a 50 MeV increase for their antiparticles, with the same 
(scaled with the number of light quarks) scenario as in i) for the baryons. Scenario i) is more suited 



for an isospin-symmetric fireball produced in high-energy collisions and was used in []15|], while 
scenario ii) may be realized at low energies. In both scenarios, the masses of the Dg mesons and 
of the charmonia are the vacuum masses. We also note that if one leaves all D-meson masses un- 
changed but allows their widths to increase, the resulting yields will increase by 11% (2.7%) for 
a width of 100 MeV (50 MeV). If the in-medium widths exhibit tails towards low masses, as has 
been suggested by |p^], to first order the effect on thermal densities is quantitatively comparable 
with that from a decrease in the pole mass. 
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Figure 6: Energy dependence of the yield of charmed hadrons relative to the charm quark pair yield for 
two scenarios of the mass change (left panel for scenario i), right panel for scenario ii), see text). For the D 
mesons, the full and open symbols are for particles and antiparticles, respectively. Note the factors 10 and 
100 for the and i/' mesons, respectively. 

The results for the two cases are presented in Fig. ^ as yields relative to the number of initially- 
produced cc pairs. As a result of the redistribution of the charm quarks over the various species, the 
relative yields of charmed hadrons may change. For example, in scenario i) the ratios of D-mesons 
are all close to those computed for vacuum masses (Fig. |5|), while for scenario ii) the changes in the 
relative abundances of the D and D mesons are obvious. In both cases the A^/D ratio is increased. 

As a result of the asymmetry in the mass shifts for particles and antiparticles assumed in sce- 
nario ii), coupled with the charm neutrality condition, the production yields of and D7 mesons 
are very different compared to vacuum masses. Overall, however, charm conservation leads to 

charmed hadrons by fixed amounts. Reducing, for example, the fight quark masses by 50 MeV wifi lower D-meson 
masses by 50 MeV and the Ac(Sc) mass by 100 (50) MeV. 
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rather small changes in the total yields. We emphasize that, although the charm conservation 
equation is strictly correct only for the total cross section we expect within the framework of the 
statistical hadronization model, also little influence due to medium effects on distributions in ra- 
pidity and transverse momentum. This is due to the fact that the crucial input into our model is 
AN^^^" /dy and there is no substantial D-meson rescattering after formation at the phase boundary. 
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Figure 7: Energy dependence of the relative change in the production yield of open charm hadrons and of 
y/i//^ meson considering different scenarios for in-medium mass modifications (see text). 

In Fig. |7| we demonstrate that the total open charm yield (sum over all charmed hadrons) 
exhibits essentially no change if one considers mass shifts, while the effect is large on charmonia. 



This is to be expected from eq. ^ : as the masses of open charm mesons and baryons are reduced, 
the charm fugacity gc is changed accordingly to conserve charm. Consequently, since the open 
charm yields vary linearly with g^ , one expects little change with medium effects in this case. In 
contrast, the yields of charmonia vary strongly, since they are proportional to g^. To demonstrate 
this we plot, in Fig. ^, the relative change of the yields with in-medium masses compared to the 
case of vacuum masses. For this comparison, we have added a third case, namely considering 
that the mass change of charmed baryons is the same as for the mesons. Because of total charm 
conservation, with lowering of their masses the open charm hadrons eat away some of the charm 
quarks of the charmonia but, since the open charm hadrons are much more abundant, their own 
yield will hardly change. 

Note that the reduction of the yield in our model is quite different from that assumed in 



0, 38, 14, 13, 16], where a reduction in D-meson masses leads to the opening up of the decay of 



and Xc into DD and subsequently to a smaller J/i/a yield from feed-down from i^' and Xc- In all 
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the previous work the in-medium masses are considered in a hadronic stage, while our model is a 
pure QGP model, with in-medium mass modifications considered at the phase boundary. 

5. Conclusions 

We have investigated charmonium production in the statistical hadronization model at lower 
energies. An interesting result is that the yield of charmed baryons (Ac, Sc) relative to the total 
cc yield increases strongly with decreasing energy. Below ^5^^= 10 GeV, the relative yield of A^ 
exceeds that of any D meson except Do, implying that an investigation of open charm production at 
low energies needs to include careful measurements of charmed baryons, a difficult experimental 
task. The charmonium/open charm yield rises only slowly from energies near threshold to reach 
~1% at LHC energy. Note that this ratio depends on the magnitude of the charm cross section, fur- 
ther underlining the importance to measure this quantity with precision. We have also investigated 
the effect of possible medium modifications of the masses of charmed hadrons. Because of a sepa- 
ration of time scales for charm quark and charmed hadron production, the overall charmed meson 
and baryon cross section is very little affected by in-medium mass changes, if charm conservation 
is taken into account. Measurable effects are predicted for the yields of charmonia. These effects 
are visible at all beam energies and are more pronounced towards threshold. 
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